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Abstract 

The paper studies the behavior of the trajectories of fluid particles in a compress- 
ible generalization of the Kraichnan ensemble of turbulent velocities. We show 
that, depending on the degree of compressibility, the trajectories either explo- 
sively separate or implosively collapse. The two behaviors are shown to result in 
drastically different statistical properties of scalar quantities passively advected 
by the flow. At weak compressibility, the explosive separation of trajectories 
induces a familiar direct cascade of the energy of a scalar tracer with a short- 
distance intermittency and dissipative anomaly. At strong compressibility, the 
implosive collapse of trajectories leads to an inverse cascade of the tracer energy 
with suppressed intermittency and with the energy evacuated by large scale fric- 
tion. A scalar density whose advection preserves mass exhibits in the two regimes 
opposite cascades of the total mass squared. We expect that the explosive sep- 
aration and collapse of Lagrangian trajectories occur also in more realistic high 
Reynolds number velocity ensembles and that the two phenomena play a crucial 
role in fully developed turbulence. 

PACS: 47.27 - Turbulence, fluid dynamics 



1 Introduction 

One of the main characteristic features of the high Reynolds number turbulent flows is a 
cascade-like transfer of the energy injected by an external source. In three dimensional 
flows, the injected energy is transferred to shorter and shorter scales and is eventually 
dissipated by the viscous friction. This direct cascade is in the first approximation 
described by the Kolmogorov 1941 scaling theory Jl] but the observed departures from 
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scaling (intermittency) remain to be explained from the first principles. As discovered 
by R.H. Kraichnan in ||, in two dimensions, the injected energy is transferred to 
longer and longer distances in an inverse cascade whereas this is the enstrophy that is 
transferred to shorter and shorter scales. Experiments || and numerical simulations 
H, H suggest the absence of intermittency in the inverse 2-dimensional cascade. In 
the present paper, we shall put forward arguments indicating that the occurrence and 
the properties of direct and inverse cascades of conserved quantities in hydrodynamical 
flows are related to different typical behaviors of fluid particle trajectories. 

Let us start by drawing some simple analogies between fluid dynamics and the theory 
of dynamical systems which studies solutions of the ordinary differential equations 

Let x s (t) denote the solution of Eq. (|1 . 1| ) passing at time s by point y. In dynamical 
systems, where the attention is concentrated on regular functions X, one encounters 
different types of behavior of solutions^ 

1) . integrable motions (more common in Hamiltonian systems), where the nearby 
trajectories stay close together forever: 

-*..»(*)! ~ 0(\ yi -y 2 \), (1.2) 

2) . chaotic motions where the distance between the nearby trajectories grows 
exponentially, signaling a sensitive dependence on the initial conditions: 

KJt) - x s> Jt)\ ~ 0(e x ^\ yi -y 2 \), (1.3) 

with the Lyapunov exponent A > 0, 

3) . last but, by no means, least, dissipative motions where 

K n (t) - x Si Jt)\ ~ 0(e x ^\ yi -y 2 \), (1.4) 

with A < 0. 

Various of these types of motions may appear in the same systems. 

Analogies between dynamical systems and hydrodynamical evolution equations, for 
example the Navier Stokes ones, are often drawn by viewing the Eulerian evolution 
of velocities as a dynamical system in infinite dimensions, see @. One has, however, 
a more direct (although not unrelated) analogy between Eq. (|1.1|) and the ordinary 
differential equation 

Tt = v(t ' x) (L5) 

for the Lagrangian trajectories of fluid particles in a given velocity field v(£, x). As 
before, we shall denote by x sr (£) the solution passing by r at time s. Clearly, the 

1 the following is not a statement about the genericity of the listed behaviors 
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system ( |1.5|) is time-dependent and the velocity field is itself a dynamical variable. 
Nevertheless, one may ask questions about the behavior of solutions of Eq. ( |1.5| ) for 
"typical" velocities. On the phenomenological level, such behavior seems to be rather 
robust and to depend on few characteristics of the velocity fields. One of them is the 
Reynolds number Re = L ^ A J- V ^ ; where |A L v| is the (typical) velocity difference over the 
distance L of the order of the size of the system and v is the kinematic viscosity. Another 
important characteristic of velocity fields is the degree of compressibility measured, for 
example, by the ratio of mean values of (J2 V a t; Q ) 2 = (V • v) 2 and ^(Vc^) 2 = (Vv) 2 . 

a a ,f3 

Reynolds numbers ranging up to O(10 2 ) are the realm of laminar flows and the 



onset of turbulence. Velocity fields in ( |1.5|) are thus regular in space and the behaviors 
(1) to (3) are observed for Lagrangian trajectories. They seem to have limited bearing 
on the character of the Eulerian evolution of velocities, see Chapter 8 of [0]. This is 
a natural domain of applications of the theory of dynamical systems to both Eulerian 
and Lagrangian evolutions. When the Reynolds number is increased, however, fully 
developed turbulent flows are produced in which the behavior of trajectory separations 
becomes more dramatic. For incompressible flows, for example, we claim (see also 
H [| [TDj] ) that the regime of fully developed turbulence is characterized by the 

2'). explosive separation of trajectories: 

l x s . ri C0 — x s ,r 2 WI becomes 0(1) in finite time. (1.6) 

More precisely, the time of separation of trajectories to an 0(1) distance is bounded 
when y 2 approaches yi, provided that the initial separation \yi — y 2 1 stays in the inertial 
range where the viscous effects may be neglected. 

Since the inertial range extends down to zero distance when Re — > oo, the fast 
separation of trajectories has a drastic consequence in this limit: the very concept of 
individual Lagrangian trajectories breaks down. Indeed, at Re = oo, infinitesimally 
close trajectories take finite time to separate^ and, as a result, there are many tra- 
jectories (in fact, a continuum) satisfying a given initial condition. It should be still 
possible, however, to give a statistical description of such ensembles of trajectories in 
a fixed velocity field. Unlike for intermediate Reynolds numbers, there seems to be 
a strong relation between the behavior of the Lagrangian trajectories and the basic 
hydrodynamic properties of developed turbulent flows: we expect the appearance of 
non-unique trajectories for Re — > oo to be responsible for the dissipative anomaly, the 
direct energy cascade, the dissipation of higher conserved quantities and the pertinence 
of weak solutions of hydrodynamical equations at Re = oo. 

The breakdown of the Newton-Leibniz paradigm based on the uniqueness of solu- 
tions of the initial value problem for ordinary differential equations is made mathe- 
matically possible by the loss of small scale smoothness of turbulent velocities when 
Re — > oo. At Re = oo, the typical velocities are expected to be only Holder continuous 
in the space variables: 

|v(t,x)-v(t,x')| 2 ~ |x-xf , (1.7) 



2 in contrast to their behavior in the chaotic regime, see Sect. 2.2 below 



3 



with the Holder exponent | close to the Kolmogorov value | The uniqueness of 
solutions of the initial value problem for Eq. ( |1.5| ) requires, on the other hand, the 
Lipschitz continuity of v(t, x) in x, i.e. the behavior (|1.7| ) with £ = 2. It should be 
stressed that for large but finite Re, the chaotic behavior (2) of trajectories may still 
persist for short separations of the order of the dissipative scale (where the viscosity 
makes the velocities smooth) and the behavior (2') is observed only on distances longer 
than that. However, it is the latter which seems responsible for much of the observed 
physics of fully developed turbulence and, thus, setting Re = oo seems to be the right 
idealization in this regime. 

For general velocity fields, one should expect that the poor spatial regularity of 
velocities v might lead to two opposite effects. On one hand, the trajectories may branch 
at every time and coinciding particles would split in a finite time as in (2'). Solving 
discretized versions of Eq. ( |1.5|) randomly picks a branch of the solution and generates 
some sort of a random walk whose average reproduces the trajectory statistics. This is 
the effect previously remarked || ^, [Kj in the studies of the incompressible Kraichnan 
model fTI|| . It should be dominant for incompressible or weakly compressible flows. On 
the other hand, the trajectories may tend to be trapped together. The most direct 
way to highlight this phenomenon is to consider strongly compressible velocity fields 
which are well known for depleting transport (see ||12|| ). An instance is provided by the 
one-dimensional equation ^| = fi(x), for (3{x) a Brownian motion in x, whose solutions 
are trapped in finite time at the zeros of (3 on the right (left) ends of the intervals where 
(3 > {p < 0) 



What may then become typical in strongly compressible velocities is, instead of the 
explosion (2'), the 

3'). implosive collapse of trajectories: 

l x s . ri CO — x s ,r 2 becomes equal to zero in finite time, (1.8) 

with the time of collapse depending on the initial positions. This type of behavior 
should lead to a domination at infinite Re and strong compressibility of shock-wave- 
type solutions of hydrodynamical equations, as in the 1-dimensional Burgers problem 



14 1 - Again, strictly speaking, we should expect behavior (3') only for Re = oo whereas 
for finite Re, on distances smaller than the dissipative scale, the approach of typical 
trajectories should become exponential with a negative Lyapunov exponent. In simple 
Gaussian ensembles of smooth compressible velocities the latter behavior and its con- 
sequences for the direction of the cascade of a conserved quantity have been discovered 



and extensively discussed in |L5| and |L6|. 



It is the main purpose of the present paper to provide some support for the above, 
largely conjectural, statements about typical behaviors of fluid-particle trajectories at 
high Reynolds numbers and about the impact of these behaviors on physics of the fully 
turbulent hydrodynamical flows. We study only simple synthetic random ensembles 
of velocities showing Holder continuity in spatial variables. Although this is certainly 
insufficient to make firm general statements, it shows, however, that the behaviors 
(2') and (3') are indeed possible and strongly affect hydrodynamical properties. In 
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realistic flows, both behaviors might coexist. For the ensemble of flows considered here, 
they occur alternatively, leading to two different phases at large (infinite) Reynolds 
numbers, depending on the degree of compressibility and the space dimension. The 
occurrence of the collapse (3') is reflected in the suppression of the short-scale dissipation 
and the inverse cascade of certain conserved quantities. The absence of dissipative 
anomaly permits an analytical understanding of the dynamics and to show that the 
inverse cascade is self-similar. This strengthens the conjectures on the general lack 
of intermittency for inverse cascades [^|, f|, [5], |T7[]. A caveat comes however from the 
consideration of friction effects, indicating that the role of infrared cutoffs might be 
subtle and anomalies might reappear in terms of them. 

Synthetic ensembles of velocities are often used to study the problems of advection 
in hydrodynamical flows. These problems become then simpler to understand than 
the advection in the Navier-Stokes flows, but might still render well some physics of 
the latter, especially in the case of passive advection when the advected quantities do 
not modify the flows in an important way. The simplest of the problems concern the 
advection of scalar quantities. There are two types of such quantities that one may 
consider. The first one, which we shall call a tracer and shall denote 6(t,r), undergoes 
the evolution governed by the equation 

d t 9 + v- WO- kV 2 6 = f, (1.9) 

where k denotes the diffusivity and f(t, r) describes an external forcing (a source). The 
second scalar quantity is of a density type, e.g. the density of a pollutant, and we shall 
denote it by p(t,r). Its evolution equation is 

d tP + V-(vp)- K V 2 p = f. (1.10) 

For / = it has a form of the continuity equation so that, without the source, the evolu- 
tion of p preserves the total mass / p(t, r)dr. Both equations coincide for incompressible 
velocities but are different if the velocities are compressible. 

A lot of attention has been attracted recently by a theoretical and numerical study 
of a model of the passive scalar advection introduced by R. H. Kraichnan in 1968 [|1^ . 
The essence of the Kraichnan model is that it considers a synthetic Gaussian ensemble 
of velocities decorrelated in time. The ensemble may be defined by specifying the 1- 
point and the 2-point functions of velocities. Following Kraichnan, we assume that the 
mean velocity (v(t, x)) vanishes and that 

(v a (t, r) v^t', r')> = 5(t - [df - d^{v - r')] (1.11) 

with constant d$ and with d af3 (r) proportional to at short distances. The latter 
property mimics the scaling behavior of the equal-time velocity correlators of realistic 
turbulent flows in the Re — > oo limit. It leads to the behavior ( |1 . 7| ) for the typical 
realizations of v. The parameter £ is taken between and 2 so that the typical velocities 
of the ensemble are not Lipschitz continuous. 

We shall study the behavior of Lagrangian trajectories in the velocity fields of a 
compressible version of the Kraichnan ensemble and the effect of that behavior on 
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the advection of the scalars. The time decorrelation of the velocity ensemble is not 
a very physical assumption. It makes, however, an analytic study of the model much 
easier. We expect the temporal behavior of the velocities to have less bearing on the 
behavior of Lagrangian trajectories than the spatial one, but this might be the weakest 
point of our arguments. Another related weak point is that the Kraichnan ensemble is 
time-reversal invariant whereas realistic velocity ensembles are not, so that the typical 
behaviors of the forward- and backward-in-time solutions of Eq. ( |1 . 5|) for the Lagrangian 
trajectories may be different. It should be also mentioned that in our conclusions about 
the advection of scalars we let Re — * oo before sending the diffusivity k to zero, i.e. 
we work at zero Schmidt or Prandtl number -. The qualitative picture should not 
be changed, however, if Re becomes very large but - stays bounded. On the other 
hand, the situation when ~ — > oo should be better described by the £ — > 2 limit of the 
Kraichnan model where the velocities become smooth. 

The paper is organized as follows. In Sect. 2 we discuss the statistics of the La- 
grangian trajectories in the Kraichnan ensemble and discover two different phases with 
typical behaviors (2') and (3'), occurring, respectively, in the case of weak and strong 
compressibility. In Sects. 3 and 4 we discuss the advection of a scalar tracer in both 
phases and show that it exhibits cascades of the mean tracer energy. In the weakly 
compressible phase, the cascade is direct (i.e. towards short distances) and it is charac- 
terized by an intermittent behavior of tracer correlations, signaled by their anomalous 
scaling at short distances. On the other hand, in the strongly compressible phase, 
the tracer energy cascade inverts its direction. In the latter case, we compute exactly 
the probability distribution functions of the tracer differences over long distances and 
show that, although non-Gaussian, they have a scale-invariant form. This indicates 
that the inverse cascade directed towards long distances forgets the integral scale where 
the energy is injected. Conversely, the scale where friction extracts energy from the 
system is shown in Section 6 to lead to anomalous scaling of certain observables. Fi- 
nally, in Section 7, we discuss briefly the advection of a scalar density. Here, in the 
weakly compressible phase, we find a cascade of the mean mass squared towards long 
distances and, on short distances, the scaling of correlation functions in agreement with 



the predictions of fll9fl . The strongly compressible shock-wave phase, however, exhibits 
a drastically different behavior with the inversion of the direction of the cascade of 
mean mass squared towards short distances. In Conclusions we briefly summarize our 
results. Some of the more technical material is assembled in five Appendices. 



2 Lagrangian flow 

The assumptions of isotropy and of scaling behavior on all scales fix the functions d a @(x) 
in the velocity 2-point function (|1.11|) of the Kraichnan ensemble up to two parameters: 

d af3 {r) = [A + (d + f - 1)B] <TV + [A-B] (r"//" 2 , (2.1) 

with A = corresponding to the incompressible case where V • v = and B = to the 
purely potential one with v = V0. Positivity of the covariance requires that A, B > 0. 
It will be convenient to relabel the constants A and B by S 2 = A+ (d— 1)B and C 2 = A. 
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S 2 and C 2 are proportional to, respectively, ((Vv) 2 ) and ((V ■ v) 2 ) and they satisfy the 
inequalities S 2 > C 2 > 0. In one dimension, S 2 = C 2 > 0. The ratio 

C 2 

P=^ 0<p<l, (2.2) 
characterizes the degree of compressibility. 

The source / in the evolution equations ( |1.9| ) and (fTTH) for the scalars will be also 
taken random Gaussian, independent of velocity, with mean zero and 2-point function 

(f(t,v)f(t',v')) = 5(t-t') X (\r-v'\), (2.3) 

where \ decays on the injection scale L. 

In the absence of the forcing and diffusion terms in Eq. ( |1.9| ), the tracer 9 is carried 
by the flow: 0(t,r) = 9(s, x tr (s)), and the density p is stretched as the Jacobian 
J = <9(x t)I .(s))/<9(r) of the map r i— > x tjI .(s): p(t,r) = p(s, x t r (s)) J. Here, x t r (s) is 
the fluid (Lagrangian) trajectory obeying dx tr /ds = t>(s,x t r ) and passing through the 
point r at time t, i.e. x ir (t) = r. The flows of the scalars may be rewritten as the 
relations 

9(t,T) = J 8(r'-^ r (a))9(s,v , )dv', p(t,r) = J 5(v - x y (t)) p(s, r') dr' (2.4) 

which imply that the two flows are dual to each other: / 9(t, r) p(t, r) dr does not change 
in time. In the presence of forcing and diffusion, there are some slight modifications. 
First, the sources create the scalars along the Lagrangian trajectories. Second, the 
diffusion superposes Brownian motions upon the trajectories. One has 



9(t,v) = 9(s,x Ur (s)) + Jf(r,x t , r (T))dT (2.5) 



for the tracer and 



p(t,v) = ^x tir ( S ))^^ + //(r,x tjr (r))^p £ iT (2.6) 



for the density, where 

^ = ^xj + v^^ x t , r (t) = r, (2.7) 
as ds 

with the overbar in Eqs. ( |2.5|) and (|2.6|) denoting the average over the (i-dimensional 
Brownian motions (3(s). 

Clearly, the statistics of the scalar fields reflects the statistics of the Lagrangian tra- 
jectories or, for k > 0, of their perturbations by the Brownian motions. In particular, it 
will be important to look at the probability distribution functions (p.d.f.'s) of the dif- 
ference r' of time s positions of two Lagrangian trajectories (perturbed by independent 
Brownian motions if k > 0), given their time t positions ri and r2, 

If (r x - r 2 , r') = ( 5(v> - x Mi (a) + (a))). (2.8) 
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Note that P 2 ' s is normalized to unity with respect to r' and the equivalent expressions 

P?{v x - r 2 , j>) = J( 5(r' + r - (*)) 5(r - ^ (s)) ) dr 

= J(5(r'-^ ri _ r (s))5(-^ r2 _ r (s)))dr, (2.9) 

where the last equality uses the homogeneity of the velocities^. 

In the Kraichnan model, the p.d.f.'s P 2 ' s may be easily computed^. They are given 
by the heat kernels e~' i_s ' M 2 (r, r') of the 2 nd -order elliptic differential operators M 2 = 
— d af3 (r) V Q V/j — 2kV 2 . What this means is that the Lagrangian trajectories undergo, 
in their relative motion, an effective diffusion with the generator M 2 , i.e. with a space- 
dependent diffusion coefficient proportional to their relative distance to power £ (for 
distances large enough that the contribution of the fc-term to may be neglected). 
Note that, due to the stationarity and the time reflection symmetry of the velocity 
distribution, 

P 2 M (r,r') = Pf(r,r') (2.10) 

but that, in general, P 2 ' s (r, r') 7^ P 2 '*(r', r), except for the incompressible case where 
the operator M 2 becomes symmetric. 

2.1 Statistics of inter-trajectory distances 

For many purposes, it will be enough to keep track only of the distances between two 
Lagrangian trajectories. We shall then restrict the p.d.f.'s P 2 ' s to the isotropic sector 
by defining 

P^(r,r')= J P*' s (Ar, r') dA = J P 2 M (r, Ar') dA, (2.11) 

SO(d) SO(d) 

where dA stands for the normalized Haar measure on SO(d). P 2 ' s (r, r') is the p.d.f. of 
the time s distance r' between two Lagrangian trajectories, given their time t distance 
r. Clearly, 

P*' s (r,r') = e-^- s]M 2(r,r') (2.12) 

with the operator M 2 restricted to the isotropic sector. In the action on rotationally 
invariant functions, 

M 2 K = -Zr^ a d r r a d r - 2Kr~ d+1 d r r d - 1 d r , (2.13) 

where 

Z = S 2 + £C 2 and a = \(d-l+i)S 2 -iC 2 }Z-\ (2.14) 

3 x t ri _ r + r coincides with x t ri in the velocity field shifted in space by r 
4 the calculation goes back, essentially, to [JOI 
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The radial Laplacian which constitutes the K-term of M£ should be taken with the Neu- 
mann boundary conditions at r = since the smooth rotationally invariant functions 
on R d satisfy d r f(0) = 0. This is the term that dominates at small r and, consequently, 
we should choose the same boundary condition^] for the complete operator M 2 . The 

oo 

adjoint operator (Mg)* with respect to the L 2 scalar product ||/|| 2 = / \f(r)\ 2 dfj,d(r), 

o 

where dfid{r) = Sd-i r d ~ 1 dr with Sa-i standing for the volume of the unit sphere in 
d- dimensions, should be taken with the adjoint boundary conditions which make the 
integration by parts possible. The diagonalization of M£ (in the isotropic sector), if 
possible, would then permit to write 

P 2 M (r,r') = fe-^SfaWM^ME), (2.15) 

where 4>e and ipE stand for the eigen-functions of the operators and (M£ )*, re- 
spectively, and du(E) for the spectral measure. We could naively expect that the same 
picture remains true for k = when 

M 2 = Ml = -Zr^ a d r r a d r (2.16) 

in the rotationally invariant sector. The problem is that the principal symbol of the 
operator M 2 vanishes at r = so that the operator looses ellipticity there and more 
care is required in the treatment of the boundary condition. 



We start by a mathematical treatment of the problem whose physics we shall discuss 
er. It will be 
transformation 



2—5 

later. It will be convenient to introduce the new variable u = r~ and to perform the 



(Uf)(u) = ( 2 -l±^f u ^-h(u^) (2.17) 

mapping unitarily the space of square integrable rotationally invariant functions on R d 
to L 2 (R + , du). The transformation U, together with a conjugation by a multiplication 
operator, turns M 2 into the well known Schrodinger operator on the half-line: 

N 2 = u~ c UM 2 U- 1 u c = Z' \-dl + (2.18) 



where 



{ ^Z, 6 = |5| and c = 6+^-1. (2.19) 



N 2 becomes a positive self-adjoint operator in L 2 (R + ) if we specify appropriately the 
boundary conditions at u — 0. The theory of such boundary conditions is a piece of 



rigorous mathematics [2C]. It says that for |6| < 1 there is a one-parameter family 
of choices of such conditions, among them two leading to the operators iVj with the 
(generalized) eigen-functions 



(p%(u) = J Tb {JE/Z'u) (2.20) 



5 this corresponds to the domination of the short distances behavior of the perturbed trajectories 
by the independent Brownian motions with the diffusion constants k 



9 



(for 6 / 0) behaving at u = as 0(u^ Tb ), respectively^} We then obtain, fixing 
the spectral measure by the dimensional consideration and e.g. the action of N 2 on 
functions vF, 



-\t-s\N?, 



±fe-M E <p%{u)<p%(v!)dE. (2.21) 



Note that the flip of the sign of b exchanges N 2 and N 2 . Relating the operators N^ 1 
to M 2 T by Eq. (pTT8|) , we infer that 

oo 

e-l^l^f^r') = ^r J e"'*— IJ-\u c <pl){r) U-\u- c <p%)(r') dE 

o 

oo 



(2-0- 



• J Tb {^E/Z'r'^) r '~ d+ * + ^ dE . (2.22) 
These are explicit versions of the eigen-function expansion (|2.15|) for k = 0. 



The eigen-functions of M 2 can be read of from the above formula. They behave, 
respectively, as 0(1) and (9(r 1_a ) at r = 0. This is the first choice that corresponds to 
the k — > limit of the Neumann boundary condition for the operator M 2 . In Appendix 
A, we analyze a simpler problem, where the operator ( |2.13| ) is replaced by its k = 
version ( |2.16| ) made regular by considering it on the interval [r , oof, with the Neumann 
boundary condition at r > 0. We show that, for |fe| < 1, this is the operator M 2 that 
emerges then in the limit r \ 0. The cutting of the interval at a non-zero value has a 
similar effect as the addition of the K-term to M 2 . We should then have the relation^ 

P 2 M (r,r') = e^l^r,/) (2.23) 

in the k — > limit, as long as \b\ < 1. 

For |6| > 1, there is only one way to make N 2 into a positive self-adjoint operator. 
If b < —1, it is still the operator N 2 that survives and the relation ( [^.23j ) still holds in 
the n — > limit. For b > 1, however, i.e. for the compressibility degree p > only 

the operator N 2 survives. Its eigen-functions (p% behave as 0(u b+J ) at u = which 
corresponds to the (9(r 1_a ) behavior of the eigen-functions of M 2 . If we impose the 
Neumann boundary condition for M 2 at r = r then, as we show in Appendix A, in the 
limit r \ 0, the eigen-functions will still become proportional to the ones obtained 
from ipg, not to those corresponding to (p^ as it happens for b < 1. The same effect has 
to occur if we add and then turn off the diffusivity k. It seems then that the equality 
P 2 s (r, r') = q~\ 1 ~ s \ m 2 (r, r') has to hold in the k — > limit when b > 1. 



6 the general boundary conditions are tp(u)\ u=0 = A u 1_2 l fc l d u u' fc '~ 7 ip{u)\ u=0 with < A < oo 

7 the choices of operators Mi with the other boundary conditions would describe the trajectories of 

particles with a tendency to aggregate upon the contact and may also have applications in advection 

problems 



10 



There is, however, one catch. A direct calculation, see Eqs. (|B.1|) and ( |B.2| ) of 
Appendix B, shows that the expression ( |2.23Q is normalized to unity with respect to r', 
but that 



-l'-l M 2 + (r,r') d» d (r') = l(b, 



r 2-e 



4Z'\t-s\ ■ 



r(6)" < i, 



(2.24) 



where 7(6, x) = f y b 1 e y dy = ^- 1 + 6; —2) is the incomplete gamma-function. 

An alternative, but more instructive, way to reach the same conclusion is to observe 
that the time derivative of e~ tM ^ (r, r') brings down the adjoint of M 2 acting on the r' 
variable so that 



J e- tM ?(r,r') r'^dr' = Z J (W^r'^ 1 ^ (■ 



Z r d r 'T 



ld-l-a+£ 



(e- tM ?{r,r') 



The contribution from r' = 00 vanishes. On the other hand, e " J " 2 (r, r'l oc r 



-tM," 



(2.25) 



tM 2 (r, r') oc r' C1+z ? , with the errors suppressed by an additional 



for small r' whereas e t}i: i>- • ' > ^ >•' "~ ri s 
factor r .It follows that the contribution from r = is zero for M 2 if 1 + a — £ > 0, 
which is the same condition as 6 < 1 or p < but it is finite for M 2 . 

The lack of normalization may seem strange since when we add the diffusivity k 
and fix the Neumann boundary conditions then, by a similar argument as for M 2 
above, the normalization is assured. The solution of the paradox is that for b > 1, the 
k — > convergence of e~* 2 (r, r') to e~ tM ? (r, r') holds only for r' 7^ and the defect 
of probability concentrates at r' = 0. For k = 0, we should then add to e~* M 2 (r, r') a 
delta-function term carrying the missing probability. 



We infer this way that 
lim Po' s (r, r) = 



for p<^, 



e-l'-W(r,r') + [1 - 7(6,3^) r(6n 5(r') 



(2.26) 



for p>4. 



In the both cases, the p.d.f.'s P 2 S satisfy the evolution equation 



d t P 2 ' s (r, r') 



T M 2 P*> s (r,r') 



(2.27) 



where M 2 , given by Eq. ( 2.16|) , acts on the r- variable and the sign =p corresponds to 
t<s. They also have the composition property: 



J pt,t ^ r ^ pt ,t (j.' jr ") dfi d (r') = Pa'* (r,r") 



ift<t'< t" or t > H > t". 
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It is instructive to note the long time behavior of the averaged powers of the distance 
between the Lagrangian trajectories. As follows from Eqs. ( |B.1| ) and (|B.2| ) of Appendix 
B, for fx > 0, 



J P2 S (r,r') r /f *d[i d (r') 



\t-s 
It - s 



2-5 



2=5 



for 



,1-a 



^ d 

P < £2: 



for p>4. 



(2.28) 



2.2 Fully developed turbulence versus chaos 

The two cases p < and p > 4s correspond to two physically very different regimes 
of the Kraichnan model. Let us first notice a completely different typical behavior 
of Lagrangian trajectories in the two cases. In the regime p < j|, which includes 
the incompressible case C 2 = studied extensively before, see fl(| and the references 
therein, the p.d.f.'s P^ir^r') possess a non-singular limit^: 

r /2-£ 

It follows that, when the time t distance of the Lagrangian trajectories tends to zero, 
the probability to find a non-zero distance between the trajectories at time s ^ t stays 
equal to unity: infinitesimally close trajectories separate in finite time. This sig- 
nals the "fuzzyness" of the Lagrangian trajectories B [nj forming a stochastic Markov 
process already in a fixed typical realization of the velocity field, with the transition 
probabilities of the process propagating weak solutions of the passive scalar equation 



d t 6 + v • V0 = [^T]]. Such appearance of stochasticity at the fundamental level seems 
to be an essential characteristic of fully developed turbulence in the incompressible or 
weakly compressible fluids. It is due to the roughness of typical turbulent velocities 
which are only Holder continuous with exponent | < 1 (in the limit of infinite Rynolds 
number Re). One should stress an important difference between this type of stochas- 
ticity and the stochasticity of chaotic behaviors. In chaotic systems, the trajectories 
are uniquely determined by the initial conditions but depend sensitively on the latter. 
The nearby trajectories separate exponentially in time at the rate given by a positive 
Lyapunov exponent. The exponential separation implies, however, that infinitesimally 
close trajectories take infinite time to separate. This type of behavior is observed in 
flows with intermediate Reynold numbers but for large Reynolds numbers it occurs only 
within the very short dissipative range which disappears in the limit Re = oo. In the 
Kraichnan model, the exponential separation of trajectories characterizes the £ —>■ 2 
limit of the fuzzy regime p < |22|, ||. 

In short, fully developed turbulence and chaos, are two different things although 
both lead to stochastic behaviors. In a metaphoric sense, the difference between the 
two occurrences of stochasticity is as between that, more fundamental, in quantum 
mechanics and that in statistical mechanics imposed by an imperfect knowledge of 
microscopic states. 



a more detailed information on how this limit is attained is given by Eq. (|B.5|) in Appendix B 



12 



2.3 Shock wave regime 

scrimp 



Let us discuss now the second regime of our system with p > 4y. In that interval, 



lim e <M 2 (r, r') = and 



lim P 2 M (r,r') = 8(r'), (2.30) 



see the second of Eqs. fl2.26p. Here the uniqueness of the Lagrangian trajectories pass- 



ing at time t by a given point is preserved (in probability). However, with positive 
probability tending to one with \t — s\ — > oo, two trajectories at non-zero distance at 
time t, collapse by time s to zero distance, as signaled by the presence of the term 
proportional to S(r') in P\ ' s (r, r'). The collapse of trajectories exhibits the trapping 
effect of compressible velocities. A similar behavior is known from the Burgers equation 
describing compressible velocities whose Lagrangian trajectories are trapped by shocks 
and then move along with them. The trapping effect is also signaled by the decrease 
with time of the averages of low powers of the distance between trajectories (< 1 — a), 
see Eq. Q2.28p , Note, however, that the averages of higher powers still increase with 



time signaling the presence of large deviations from the typical behavior. 

Due to the inequalities < p < 1, the second regime, characterized by the collapse 
of trajectories, is present only if £ 2 > d, i.e. for d < 4. Its limiting case with £ = 2 and 
d < 4 was first discovered and extensively discussed in ]nj and flEf . It appears when 



the largest Lyapunov exponent of (spatially) smooth velocity fields becomes negative. 



3 Advection of a tracer: direct versus inverse cas- 
cade 

3.1 Free decay 

Let us study now the time t correlation functions of the scalar 9 whose evolution is given 
by Eq. ( |1.9|) . Assume first that we are given a homogeneous and isotropic distribution 
of 9 at time zero and we follow its free decay at later times. From Eqs. fl2.4|) and (|2.9| ) 
we infer that 

oo 

F*(fi,r) = (9(t,r)9(t,0)) = J P*'\r,r') F|(0, r') d^ d (r'). (3.1) 

o 

In particular, to calculate the mean "energy" density e g (t) = (?9(t,r) 2 ) = ^F|(t,0), 
the separation r should be taken equal to zero. For p < |r, the limit P^^O, r') is 
a regular positive function and it stays such even for k = 0, see Eq. ( |2.29| ). Since 
F|(0, r') < F|(0,0) as a Fourier transform of a positive measure, it follows that the 
total energy diminishes with time: e g (t) < e e (0). On the other hand, for p > ^, 

P 2 *'°(0,r') = 5(r'), see Eq. ( ^30) ), and the total energy is conserved: e e (t) = e e (0). 
The loss of energy in the regime p < is not due to compressibility^, but to the 

9 in temporally decorrelated velocity fields, the mean energy e g is conserved also in compressible 
flows, in the absence of forcing and diffusion 
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non-uniqueness of the Lagrangian trajectories responsible for the persistence of the 
short-distance dissipation in the k — ► limit. As is well known this dissipative anomaly 
accompanies the direct cascade of energy towards shorter and shorter scales in the 
(nearly) incompressible flows. On the other hand, in the strongly compressible regime 
p > the scalar 9 is transported along unique trajectories and its energy is conserved 
in mean. The short distance dissipative effects disappear in the limit k — > 0: there is 
no dissipative anomaly and no direct cascade of energy. As we shall see, the energy 
injected by the source of 9 is transferred instead to longer and longer scales in an inverse 
cascade process. 



3.2 Forced state for weak compressibility 

The direction of the energy cascade may be best observed if we keep injecting the energy 
into the system at a constant rate. Let us then consider the advection of the tracer in 
the presence of stationary forcing. From Eqs. (|2.5| ) and (|2.3| ), assuming that 9 vanishes 
at time zero, we obtain, 



oo 

t.s, 



m,r) = (Jf(s,^ r (s))dsff(s',^ (s'))ds') = ds P}' s (r,r')x(r')dn d (r'), (3.2) 

J J 



which is a solution of the evolution equation 

d t F e 2 = -M*Fl + X (3-3) 

with the operator M 2 given by Eq. (|2.13|) . 

When p < d ~2^ (i.e. for a > 1 or b < 0), which implies that we are in the weakly 
compressible phase with p < t|, and for k = 0, 



_ t oo oo 

F°(t,r) = ids! e~ sM * (r,r') X (r') d^(r') f (M.^)' (r,r') X (r') dfi d (r') 



= 7 r' 1 ^ X(r) dr> + ^ r 1 ^ / r' a ^ X {r>) dr' = F e 2 {r) , (3.4) 

see Eq. ( |C.5|) of Appendix C. Thus for p < d ~ 2 2 g +g , when t — > oo, the 2-point function of 9 
attains a stationary limit F® (r) with a finite mean energy density e e = (\9 2 ) = JF|(0). 
The corresponding stationary 2-point structure function is 

S e 2 (r) = ((9(r)-9(0)) 2 ) = 2(F 2 e (0) - F e 2 {r)) = flC a dCIr' a ^ X (r')dr' 

o o 

(2-OU + a-OZ ^ fOT r SmaU ' 

(a _ 2 1)z / r jl ~^ x{ r ') dr' — (a _ 2 i )z r-^- 1 ) J r /a ~? x (r') dr' for r large. 
o o 

Thus S 2 (r) exhibits a normal scaling at r much smaller than the injection scale L 
whereas at r > L the approach to the asymptotic value 2(9 2 ) is controlled by the 
scaling zero mode r 1_a of the operator M 2 . In Appendix D, we give the explicit form 
of the stationary 2-point function F 2 in the presence of positive diffusivity k. 
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3.3 Dissipative anomaly 

Let us recall how the dissipative anomaly manifests itself in this regime. The stationary 
2-point function of the tracer solves the stationary version of Eq. ( |3.3|) . When we let 
in the latter r — > for positive k, only the contribution of the dissipative term in 
M% survives and we obtain the energy balance equation e e = k ((V0) 2 ) = ix(0)> i- e - 
the equality of the mean dissipation rate e e and the mean energy injection rate ix(0)- 
Taking first k — > and r — > next, instead, we obtain the analytic expression of the 
dissipative anomaly: 

lim e e = \ lim M 2 lim F 9 2 {r) = i x (0). (3.6) 

Thus, in spite of the explicit factor k in its definition, the mean dissipation rate does 
not vanish in the limit k — > 0, which explains the name: anomaly. 

For ^ti < p < (1 e for < 6 < 1, see Eq. (|CT|)), the 2-point function F°(t, r), 
still given for n = by the left hand side of the relation ( |3.4|) , diverges with time as t b . 
More exactly, as we show in Appendix C, it is the expression 

*Z(t,r) - (1 Jgg;_ t) t» 7r' a -' X (r')dr> (3.7) 

that tends to the right hand side of the relation (|3.4j). Finally, for p = d ~ 2 2 g +e , there is a 
constant contribution to F|(t,r) logarithmically divergent in time. For d ~ 2 2 s +s < p < 
the system still dissipates energy at short distances with the rate e g that becomes equal 
to the injection rate asymptotically in time, but it also builds up the energy e g (t) in the 
constant mode with the rate decreasing as t~^~ b \ Note that in spite of the divergence 
of the 2-point correlation function, the 2-point structure function of the tracer still 
converges as t — > 00 to a stationary form given by Eq. ( |3.5|) . Now, however, S^r) is 
dominated for large r by the growing zero mode oc r l ~ a of M2. 

3.4 Forced state for strong compressibility 

Let us discuss now what happens under steady forcing in the strongly compressible 
regime p > ■£ (i.e. for 1 + a — £ < or 6 > 1). Here the 2-point function ( |3.2|) , which 
still evolves according to Eq. (|3.3|), is for k = given by the relation 

t 00 _i_ 
F°(t, r) = fdsf e~ sM * (r, r') X (r') dfx d (r') 


+ x(0)Jil- 1 (b,g^)T(b)- 1 ]ds. (3.8) 

When t — > 00, the first term on the right hand side tends to 

00 

/(M 2 + )- 1 (r,r')x(r') dfM d (r') 


= (I^jz / r ' 1_€ *( r ') dr> + jd^z ^ J r '^ dr ' ( 3 - 9 ) 
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/ r' 1 + a -Z x '(r')dr' 
(2-C)(l + a-C)Z r ~~ (l-a)(l + a-OZ r 101 r SIlla11 5 

(3.10) 



00 i-f 

_L_ J r x ( r ') ^ r f or r l arge; 


where the asymptotic expressions hold for 1 + a — £ < 0, i.e. for p > On the other 
hand 

X(0) } (b r»-f )d X(0) 2-g f311 N 

r(6) 4Z'(t-,)J flS ^ (2-0(1+0-0^ ' 1 j 

except for p = j| when it diverges as x< ^\ r z , Int. Hence, for p > |j, the quantity 
-F|(t, r) — x(0)t converges when t — > oo and the limit is proportional to the zero mode 
r 1_a of M 2 for small r (up to 0(r A ~^) terms). As we see, the energy injected into the 
system by the external source is accumulating in the constant mode with the constant 
rate equal to the injection rate hx(fy- The dissipative anomaly is absent in this phase. 
Indeed, 

e s = ~ lim lim M 2 F|(t,r) = (3.12) 

and the same is true at finite times since F$(t, r) becomes proportional to the zero 
modes of M 2 at short distances. These are clear signals of the inverse cascade of energy 
towards large distances, identified already in the £ — ■> 2 limit of the p > -A regime in 



n mi- 

The 2-point structure function 

S e 2 (t,r) = 2(F|(t,0) -F 2 e (t,r)) = 2 X (0)t - 2F°(t,r), (3.13) 

which satisfies the evolution equation 

dtS* = -M 2 S 6 2 + 2( X (0) -x), (3-14) 

reaches for p > the stationary limit whereas it diverges logarithmically in time for 
p = |y. Note that it is now at large r that S%(r) scales normally oc r 2 ~^ and at small r 
that it becomes proportional to the zero mode r x ~ a of M 2 . 



4 Intermittency of the direct cascade 

The higher correlation functions of the convected scalars involve simultaneous statistics 
of several Lagrangian trajectories. To probe deeper into the statistical properties of the 
trajectories, it is convenient to consider the joint p.d.f.'s P^ s (ti, . . . ,r N ; r[, . . . ,r' N ) of 
the time s differences of the positions r' 1; . . . of N Lagrangian trajectories passing at 
time t through points ri, . . . , r^. In the notation of Section 2, 

r N 

P t /(r l ,...,r N ;v' l ,...,v' N ) = / ( J] S(r' n - x t>r „( S ) + r) ) dr. (4.1) 

J n=l 
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Clearly, the functions -P/v S (r;rQ are translation- invariant separately in the variables 
r (ri, . . . , r N ) and in = (r' 1; . . . , r' ). In the Kraichnan model, the p.d.f.'s P t,s are 



again given by heat kernels of degree two differential operators |23 

i^ferO = e-^Hr-l), (4.2) 

where the operators 

K = E d^(r n -r m )V r ,V ri - « £ V r 2 „ (4.3) 

l<n<m<V l<n<iV 

should be restricted to the translation-invariant sector, which enforces the separate 
translation-invariance of their heat kernels. The relations -P^' s (r;rQ = P*'*(r; rQ gen- 
eralize Eq. ( 2.10|) . As for N = 2, P*' s (r; = P*'*(r^; r) only in the incompressible 
case. 

Strong with the lesson we learned for two trajectories, we expect completely different 
behavior of the p.d.f 's P^ s {r; rQ for r n 's close to each other in the two phases, resulting 
in different short-distance statistics of convected quantities. Let us start by discussing 
the weakly compressible case p < j^. Here we have little to add to the incompressible 
story, see e.g. 0, |IO|j . We expect the limit lim P^ s (r; rQ = P^ s (0; rQ to exist and to 

be a continuous function (except, possibly at = 0) decaying with \t — s\ and at large 
distances, just as for P 2 ' s , see Eq. Q2.29| ). More exactly, we expect || an asymptotic 
expansion generalizing the expansion (|R3l) of Appendix B for P*' s : 



P'/i^X) = E * ai+{2 ~° j <ki®M\t-*\,& (4-4) 



N 

i 

3=0,1,.. 



for A small, where 4>i,o are scaling zero modes of the operator = M N with scaling 
dimensions oi > and <pi )P are "slow modes", of scaling dimension <7j+(2— satisfying 
the descent equations M N (pij = (pij-i- The constant zero mode 0o,o = 1 (corresponding 
t° V'o.o = -Pjv(Q; ■ )) gives the main contribution for small A, but drops out if we consider 
combinations of P t,s (Ar; rQ with different configurations r which eliminate the terms 
that do not depend on all (differences) of r„'s. Such combinations are dominated by 
the zero modes depending on all r n 's. For small £, there is one such zero mode <pi 0j o for 



each even N. A perturbative calculation of its scaling dimension done as in |26|], where 
the incompressible case was treated, gives 



N 

In the absence of forcing, the iV-point correlation functions F e (t, r ) = ( n @(t, r n ) ) 

n=l 

of the tracer are propagated by the p.d.f.'s P*' s : 

KM = JPi s (rr,l)F N (s,T)d'r! (4.6) 
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where dV = dr' 2 - ■ ■ dr' , compare to Eq. (|3.1|) . In the presence of forcing, F e obey 



recursive evolution equations [24, 23 1. If F vanish at time zero then the odd correlation 



functions vanish at all times and the even ones may be computed iteratively: 

t 

F e N (t,r) = fdsf Pj'ferQ £ Fl^v^.^.^r'J X (\r' n - r'J) dV. (4.7) 



We expect that for small £ or weak compressibility, F e (t, r) tend at large times to the 
stationary correlation functions F e (r) whose parts depending on all r n 's are dominated 
at short distances by the zero modes of M N . In particular, this scenario leads to the 
anomalous scaling of the iV-point structure functions S 9 N (r) = {(6(r) — 6(0)) N ) which 
pick up the contributions to F e N dependent on all r n 's. Naively, one could expect that 



S e (r) scale for small r with powers y(2— £), i.e. y times the 2-point function exponent. 
Instead, they scale with smaller exponents, which signals the small scale intermittency: 



S° N {r) oc r ( 2 ~0*+^ (4.8) 
with the anomalous (part of the) exponent A e N given for small £ by 

K = - N{N ;Z) 2P) * + °^ 2 )' ( 4 - 9 ) 

see Eq. (|4.5| ). We infer that the direct cascade is intermittent. 



5 Absence of intermittency in the inverse cascade 
5.1 Higher structure functions of the tracer 

For p > i.e. in the strongly compressible phase, we expect a completely different 
behavior of the p.d.f.'s P^ S {L', iL) when the points r n become close to each other. The 
(differences of) Lagrangian trajectories in a fixed realizations of the velocity field are 
uniquely determined in this phase if we specify their time t positions. The p.d.f.'s for 
iV trajectories should then reduce to those of iV — 1 trajectories if we let the time t 
positions of two of them coincide: 

hm PJ-CeeO = P^^,...,^;^,...,^)^-^). (5.1) 

1 JV 1 N-1 

N 

Applying this relation N times, we infer, that P' ,s (r, . . . , r; rQ = T7 <K r i — r n) • Since 

n=2 

the p.d.f.'s P t,s propagate the A^-point functions of the tracer in the free decay, see 
Eq. (|4.6|) , it follows that, in the strongly compressible phase, such a decay preserves 
all the higher mean quantities (0(t, r) N ) = F^(t, r, . . . , r). In the presence of forcing, 
however, all these quantities are pumped by the source. Indeed, Eq. (|4.7|) implies now 
that 

(0{t, rf) = ^ X (0) Ms, O^ 2 ) da, (5.2) 
z o 
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from which it follows that, for even N, (9(t, r) N ) = (N — 1)!! (x(O)t) . 



The relation ( |5.1|) permits also to calculate effectively the higher structure functions 
S e N (t,r) in the strongly compressible phase. We prove in Appendix E that for N even, 

S N (t,r) = N(N-l) } dsj Pt>\r,r') S e N 2 {s,r') ( x (0) - X (r')) d^r'). (5.3) 

oo 

Note that S e N satisfies the evolution equation 

d t S d N = -M 2 S d N + N(N - 1) S° ( X (0) - X ). (5-4) 



This is the same equation that would have been obtained directly from Eq. (|1.9|) ne- 
glecting the viscous term and averaging with respect to the Gaussian fields v and /, 



e.g. by integration by parts ||25|| . The situation should be contrasted with that in the 
weakly compressible case where the evolution equations for the structure functions do 
not close due to the dissipative anomaly which adds to Eq. (|5.4|) terms that are not 



directly expressible by the structure functions [JO]], see also |26| . 



Substituting into Eq. (|5.3|) the expression (|2.26|) for P 2 ' s in the strongly compressible 
phase, we obtain 



t oo 



S 9 N (t,r) = N(N-1) JdsJe-^ M ^r,r')S d N a (a/) (x(O)-x(r')) WO- (5-5) 

oo 

The above formula implies, by induction, that S e N are positive functions (no surprise), 
growing in time. Suppose that S e N _ 2 (t,r) reaches a stationary form (r) which 
behaves proportionally to the zero mode r l ~ a for small r and which exhibits the normal 
scaling oc r^-iX 2- ?) for large r (S^ behaves this way for p > i.e. for b > 1). Then 

Jds7e-^ M ^r,r') S e N 2 (s,r') X (r') dfi d (r') (5.6) 
o o 

converges when t — > oo to a function bounded by 

oo 

/(M+)- 1 (r,r') S e N 2 (r') X W) dfx d (r'), (5.7) 
o 

which behaves as r 1_a for small r and as a constant for large r, compare to the estimate 
(|3.10|) . On the other hand, the dominant contribution to the x(0) term in Eq. (|5.5| ) is 
proportional to 



/ ds Je- sM * (r, r') r /( ^ 1)(2 ^ } d^ d {r') . {h.i 
o o 



Since, by Eq. ( p.2|) of Appendix B, / e sM i (r, r') r'^ ^ dfid{f') vanishes at s = 

o 

and behaves as ,g"2"~ 1-{ y 1-a for large s, we infer that the integral (|5.8|) stabilizes when 
t — > oo only if 

f < b = g . (5.9) 
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This condition becomes more and more stringent with increasing N. If it is not satisfied, 
then the contribution ( |5.8| ), and, consequently, S e N (t,r), diverge as t^"~ r 1_a . If it is 
satisfied, the contribution (|5.8|) reaches a limit when t — > oo which is proportional to 
r -r( 2 -€) _ then dominates for large r the stationary Appoint structure function S e (r) 
which for small r behaves as r 1_a , reproducing our inductive assumptions. 

Summarizing: The even iV-point structure functions become stationary at long 
times only if the conditions ( |5.9|) are satisfied and they exhibit then the normal scaling 
at distances much larger than the injection scale L, i.e. in the inverse energy cascade. At 
the distances much shorter than L, however, the existing stationary structure functions 
are very intermittent: they scale with the fixed power 1 — a. 



5.2 Generating function and p.d.f. of scalar differences 

It is convenient to introduce the generating function for the structure functions of the 
scalar defined by 



OO 



Z\\t,r) = { ^m,r)-e m)) = J2t^fl S L(t,r). (5.10) 

n=0 

We shall take A real. Note that the evolution equation Eq. ( |5.4j) implies that 

d t Z e = _[M 2 + A 2 ( X (0)-x)]^. (5.11) 

At time zero, Z = 1. Since M 2 + A 2 (x(0) — x) nas a similar boundary condition problem 
at r = as M 2 , one should be careful writing down the solution of the parabolic equation 
( |5.11| ). It is not difficult to see that Z e is given by a Feynman-Kac type formula: 

Z e (X,t,r) = E r (e o ), (5.12) 

where E r is the expectation value w.r.t. the Markov process r(s), s > 0, with transition 
probabilities -P 2 ' s (r, r'), starting at time zero at r. Due to the delta-function term in 
the transition probabilities, almost all realizations r(s) of the process arrive at finite 
time at r = and then do not move. Note that Z e (0,t,r) = Z e (X,t,0) = 1 and that 
Z e (X,t,r) decreases in time. Moreover, Z e (X,s,r) for s > t is bounded below by an 
expectation similar to that of Eq. ( |5.12[ ) but with the additional restriction that r(t) = 
(and hence that r(s') = for all s' > t). The latter is positive since the probability 
that r(i) = is non-zero (it even tends to one when t — > oo). Thus a non-trivial limit 
lim Z 9 (X,t,r) = Z 9 (X,r) exists. It satisfies the stationary version of Eq. ( |5.11| ): 



[M 2 + A 2 (x(0)- X )]^ = 0. (5.13) 

In particular, for r large, for which we may drop x( r ); Z e (X,r) is an eigen-function of 
the operator M 2 given by Eq. ( j2.16 ) with the eigen- value — A 2 x(0). This permits to 



find the analytic form of the generating function Z (A,r) in this regime, a rather rare 
situation in the study of models of turbulence. We have 

Z\\r) = ^ ^ X (0)/Z'\X\r^ b K b (^ X (0)/Z'\X\r^) = Z e sc (X 2 r^). (5.14) 
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The Bessel function Ki(z) decreases exponentially at infinity. We have chosen the 
normalization so that Z 9 (0,r) = 1. Since z h Ki ) {z) has an expansion around z = in 
terms of z 2n and z 2b+2n with n > 0, it is iV-times differentiate le at zero only if y < b. 
Not surprisingly, this is the same condition that we met above as the requirement for 
the existence of stationary limits of the structure functions S 9 N (t,r). 

In short: in the strongly compressible phase p > j|, the generating function 
Z 9 (X,t,r) has a stationary limit Z e (X,r) which for large distances takes the scaling 
form Z 9 c (X 2 r 2 ~^). Although Z 8 (X,r) is non-Gaussian and even not a smooth function 
of A, its normal scaling in the large r regime, responsible for the normal scaling of 
the existing structure functions, signals the suppression of intermittency in the 
inverse cascade. 

The same behavior may be seen in the Fourier transform of the generating function 
Z e (\,t,r) giving the p.d.f. of scalar differences: 

V 9 {t,$,r) = (6(#-6(t,T)+9(t,0))) = ± Je- iX »Z 9 {\,t,r)d\. (5.15) 

The t —>■ oo limit V 9 ^^) of the finite-time p.d.f. satisfies the partial-differential 
equation [M 2 — (x(0) — x) d 2 } V 9 = 0. For r large, the latter reduces to the ordinary 
differential equation 

d x [( X (0) + Z'x 2 )d x + (2b+l)Z'x]p = 0, (5.16) 
where V 9 ($, r) = r~^~ p 9 {r~^~ $) . The normalized, smooth at x = solution is 



A') - 1L X -'y [xW + ZVr 1 -*. (5.17) 

It is the Fourier transform of the generating function Z 9 C (X 2 ). Note that the condition 
y < b becomes now the condition for the existence of the iV-th moment of p 9 . The 
slow decay of p 9 (x) at infinity renders most of the moments divergent. 



6 Infrared cutoffs and the inverse cascade 

As shown in the previous section, in the strongly compressible phase, the asymptotic 
behavior of the scalar 9 is quasi-stationary: due to the excitation of larger and larger 
scales, observables might or might not reach a stationary form. It is therefore of fun- 
damental interest to analyze how the inverse cascade properties are affected by the 
presence of an infrared cutoff at the very large scales. This has also a practical impor- 
tance as such cutoffs are always present in concrete situations in one form or another 
[0, |], H l7j . The simplest modification of the dynamics that introduces an infrared 



cutoff is to add to ( |1.9| ) a friction term: 

d t 6 + v -V0 + a6 - kV 2 6 = f, (6.1) 
where a is a positive constant. We shall be interested in studying the limit when a — > 0. 



For flows smooth in space and 5-correlated in time, the case considered in [27|, the 
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advection and the friction terms have the same dimensions. For non-smooth flows, this 
is not the case and friction and advection balance at the friction scale r\ f ~ O faT 1 ^ 2- ^! 
which becomes arbitrarily large when a tends to zero. Roughly speaking, advection and 
friction dominate at scales much smaller and much larger than rj f , respectively. The 
hierarchy of scales is therefore the mirror image of the one for the direct cascade : the 
energy is injected at the integral scale L, transferred upwards by the advection term 
in the convective range of the inverse cascade and finally extracted from the system at 
very large scales. We are interested in the influence of the infrared cutoff scale on the 
inverse cascade convective range properties and we shall therefore assume that p > ^ 
throughout this section. 

Heuristically, it is a priori quite clear that the friction will make the system reach a 
stationary state. Specifically, the friction term in Eq. ( |6.1| ) is simply taken into account 
by remarking that the field 9(t, r) = exp(a t)9(t, r) satisfies the equation (|1.9|) with a 
forcing e at f(t,r). We can then carry over the Lagrangian trajectory statistics from 
the previous Sections and we just have to calculate the averages with the appropriate 
weights. For example, the recursive equation ( |4.7| ) for the iV-point function in the 
presence of forcing becomes 

t 

F e N (t,r) = fdsj e-(*-)^Pj-(r;rO £ F°_ 2 (s, v[, ... ., r^) x(|r' n -r'J) rfV. (6.2) 



Similarly, the expressions ( |5.2j ) and ( |5.3| ) for (6(t,r) N ) and S e N (t,r) are modified by 
inclusion of the factor e -( t ~ s ) Na under the time integrals. This renders them convergent 
in the limit t — > oo, in contrast with the a = case. As a result, (0(t, r) N ) and S e (t, r) 
reach when t —>■ oo the limits that are the solutions of the stationary versions of the 
evolution equations 

d t (9 N ) = &i! x (0)(^- 2 > _ Na(9 N ), (6.3) 
d t S e N = -M 2 S e N + N(N - 1) S e N _ 2 ( X (0) - x )-NaS N . (6.4) 

, . N 

We obtain then in the stationary state: (0(r) ) = (N — 1)!! (^^) T and 



A JV 
~? ~2 



v y k=l v 7 

where the variables Zk are defined as 

Zk = 2kH^ with V f = C~^)^ (6.6) 



being the friction scale. For r ^> rj , all the Bessel functions tend to zero and S 8 N reaches 

N N 

the constant asymptotic value 2 (N— 1)!! ( 2s ^r)^~ which agrees with the stationary value 

N 

of 2 T (6 N ). The expansion of for small arguments gives the following dominant 
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Figure 1: The probability distribution function of scalar differences with and without 
friction (solid and dashed lines). The specific parameters are d=l,!; = l,S 2 =C 2 = j, 
X (Q) = l ; a = 2, andr = 0.01. 



behaviors in the inverse cascade convective range L^r-C^: 



C\ r 
c 2 r 



f (2-0 



if b > f 



f (2-0 



■<*(*) 



if 6 



2 • 



2 • 



f (2-0 /^N (2-0(^-6) 



, c 3 r - (^) 
where the constants Cj are given by 



if 6<f, , 



(6.7) 



Ci 



^(0)^ m r(b-AT/2) 
V 4Z' y (AT/2)! r(6) ' 



c 2 



V 4Z' y (AT/2)! r(6) ' 



° 3 V 4Z' y (iV/2-1)! r(l+6) I 4^ fc v fc-i / 



(6.* 



The threshold 6 = ~ in Eq. (|6.7|) is the same as in the inequality (|5.9|) , discriminat- 
ing the moments that do not converge at large times in the absence of friction. The 
converging moments are not modified by the presence of friction. Conversely, those 
that were diverging are now tending to finite values but they pick an anomalous scaling 
behavior in the cutoff scale rj f . Note that the moments with ~ > b scale all with the 
scaling exponent 1 — a. 

It is interesting to look at this saturation from the point of view of the p.d.f. r) , 

defined as in the relation ([5.15 ). The equation for V e can be derived by the same 
procedure as in the previous section. Its stationary version reads 



- M 2 V e + a d<> ($V ) + (x(0) - x) dlV e = 0. 



(6.9) 
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For the scales r > L of interest to us here, \ can be neglected with respect to x(0)- 
The relevant informations on the p.d.f. V e are conveniently extracted by expanding the 
function V e in a series 



rfrr) = ^e-«E4«^ [iJmV (6 - 10) 

of the Hermite polynomials H 2 j,. The coefficients d 2 k can be obtained by plugging the 
expansion ( |6.10|) into Eq. (|6.9| ) and using well-known properties of Hermite functions. 
One obtains 

d {r) = l, d 2k (r) = k]2 b-i +2 k A K b (z k ), (6.11) 

where z k is defined in Eq. ( |6.6|) . At scales r rj f , friction dominates and the p.d.f. 
tends to a Gaussian form. In the inverse cascade convective range L ^ r ^ t] f , the 
solution ( |5.17| ) remains valid as long as $ 2 <C while the power-law tails are cut by 
an exponential fall off. The situation is exemplified in Fig. 1. The scaling behavior (|6.7f) 
of the structure functions is then easy to grasp: for y < b, the dominant contribution 
comes from the scale invariant part of the p.d.f., resulting in the normal scaling, whereas 
for ~ > b, the leading contribution comes from the region around 9 2 = with the 
tails cut out by friction. The dominant behavior can be captured by simply calculating 
the moments with the scale- invariant p.d.f. ( [5.17D cut at $ 2 - 



a 



7 Advection of a density 

The time t 2-point function of the scalar p, whose advection is governed by Eq. (|1.10D , 
may be studied similarly as for the tracer 6, see |T9| , 28]. For the free decay of the 
2-point function, we obtain 

oo 

F^t,r) = (p(t,v)p(t,0)) = Jp^(r',r)F 2 p (0,r')dfi d (r'), (7.1) 



where we have used Eqs. ( ]2.4| ) and ( p.9| ), compare to Eq. ( |3.1| ). The evolution ( |1.10| ) 



of the scalar p preserves the total mass Jp(t,r)dr. As a consequence, the mean "total 

oo 

mass squared" per unit volume m 2 (t) = J {pit, r)p(t, 0))dr = j F 2 (t, r) dpdij") does not 

o 

change in time in both phases. 

In the presence of the stationary forcing, the 2-point function of p computed with 
the use of Eqs. ( p.6|) and ( ]2.9|) becomes 

t oo 

Fg(t,r) = JdsJ Pf(r',r) X (r') dp d (r'), (7.2) 
o o 

if p = at time zero. It evolves according to the equation 

d t F 2 p = -(M«)*it + X, (7-3) 
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i.e. similarly to _F|, see Eq. (|3.3| ), but with M£ exchanged for its adjoint (M£)*, a 
signature of the duality between the two scalars noticed before. 

For p < d ~ 2 2 f +g , the 2-point function F£(t, r) attains a stationary form given by 
Eq. ( p. 2D of Appendix D for k > and reducing for k = to the expression 

F ^ r ) = J^I)z r ~ d+1+a ^ X(r') dr> + ^r-Wfr'"- 1 X {r>) dr> . (7.4) 

In particular, F^r) becomes proportional to < r - d + 1 + a -z for small r and diverges at 
r = 0, except for the incompressible case when the two scalars coincide. The small r 
behavior agrees with the result of |1| and, in one dimension, with that of [p8fl . For 
large distances r, the function F^r) is proportional to r~ d+2 ~^. In the upper interval 
d ~2^ < p < fr of the weakly compressible phase, it is 

Wt, r) - (.J, 4 ^, I* r-«+-< / r"» *(r<) 

that reaches the stationary limit still given for k = by the right hand side of Eq. ( |7.4j ). 
Thus the 2-point function F^t, r) is pumped now into the zero mode r ~ d+1+a ~£ of the 
operator M^. 

N 

The higher order correlation functions of p, F^(t,r) = ( T7 p(^ r n)), are expected 

n=l 

to converge for long times to a stationary form for sufficiently small £ and/or compress- 
ibility and to be dominated at short distances by a scaling zero mode ipo of the operator 

ML) = l -dpi Y, ln l r n-r m | + O^ 2 ). (7.5) 

l<n<m<N 

The zero mode ipo becomes equal to 1 when £ — > 0. The scaling dimension of ipo may 
be easily calculated to the first order in £ by applying the dilation generator to the left 
hand side of Eq. ( |7.5|) . It is equal to 

"Wpt +0(e) = £(2-0 + A' (7.6) 



2 



which again agrees with the result 

AP = -N + JV(1 ~ (JV ~ 1)pd) £ + Q(i 2 ) 



of JT9 and with the exact result = a — 1 — e£ obtained above. Note the singular 
behavior of ^od) at the origin, at least for small £. 

Finally, in the strongly compressible phase p > j^, where the second of the expres- 
sions ( p.26| ) has to be used for P 2 S '*, we obtain 

t OO _L 

F 2 p (t,r) = fdsf e~ sM * (r',r) x (r') dptf) 
o o 
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t OO ,0 — p 

+ 5(r) Jds J [I- 7 (6, r^.) Tib)- 1 } X {r') e^(r') . (7.7) 
oo 

For large times, F$(t,r) is pumped into the singular mode proportional to the delta- 
function with a constant rate0: 

F£(t,r) - 8{v)t J X {r')dUr') 7^^+^] r> d - a X {r>) dr> 

+ (T=W r ^ + " ? / ^ dr ' + 5 « X(0 dr> (7.8) 



(except for p = ^), compare to Eqs. (|3.9|) and (|3.11| ). Its non-singular part, however 



- j 

stabilizes and becomes proportional to r~ d+2 ~^ for small r and to r~ d+1+a ~^ for large r. 
Note the inversion of the powers as compared with the weakly compressible phase. 

The mean total mass squared of p per unit volume, m 2 , exhibits in the presence 
of forcing a position-space cascade analogous to the wavenumber-space cascade of the 
energy e e . Let us localize m 2 in space by defining its amount between the radii r and 
R as 

KrAt) = ?F 2 p (t,r')dp d (r'). (7.9) 



Integrating the evolution equation (|7.3|) for F£ with the radial measure dpd from r to R, 
and using the explicit form of = — Z r~ d+1 d r r a d r r d ~ l ~ a+ ^ , we obtain the relation 

dtm% rtR (t) = ZS d ^ 1 r' a d rl r' d - 1 - aH F 2 p (t,r') R + f X (r') d^r') (7.10) 

r r 

expressing the local balance of the total mass squared, provided that we interpret 

R 

/ x{ r ') d[id{ r ') as the injection rate of m 2 in the radii between r and R and 

r 

ZSd^dy- 1 -^ F 2 p (r,t) = $(r) 
as the flux of m 2 into the radii < r. In the weakly compressible phase p < 4- and in the 

r 

stationary state, $(r) = — / X {r') dpd{ r ') so that the flux is constant for r much larger 

o 

than the injection scale L and it is directed towards bigger radii. On the other hand, in 

oo 

the strongly compressible phase p > 4-, one has the equality <3>(r) = / X {r') dfj,d(r') so 

r 

that the flux is directed towards smaller distances. It eventually feeds into the singular 
mode all of m 2 injected by the source. As we see, the two phases differ also by the 
direction of the cascade of the total mass squared of p. 

8 Conclusions 

We have studied the Gaussian ensemble of compressible d- dimensional fluid velocities 
decorrelated in time and with spatial behavior characterized by the fractional Holder 



10 the pumping disappears if /o° x( r )^/ i rf( r ) = 0; which is the case considered in [£8|, 19 , but even 
then F£ picks up a singular contribution in the limit K — > 
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Figure 2: The second-order structure function S^ir) vs r for £ = 1.8 and d = 3 in the 
three different regimes b = — 2 (dot-dashed line), b = 0.5 (dashed line) and 6 = 2 (solid 
line). 

exponent |. We have shown that the Lagrangian trajectories of fluid particles in such 
an ensemble exhibit very different behavior, depending on the degree of compressibility 
p of the velocity field. For p < i.e. b defined in (|2.19| ) smaller than unity, the 
infinitesimally close trajectories separate in finite time, implying that the dissipation 
remains finite in the limit when the molecular diffusivity k — + and that the energy is 
transferred towards small scales in a direct cascade process. The constancy of the flux 
at small scales leads to a normal scaling behavior r 2 ~^ of the second order structure 
function 5f(r) for r C I (the typical scale where the energy is injected). For b negative 
(which includes the incompressible case), as the system evolves, the dissipation rate 
tends to the injection rate rapidly enough to ensure that the energy (6 2 ) remains finite. 
The non-constant zero mode r^ 2 ~^ b controls the decay of (r) to its finite asymptotic 
value 2(6 2 ) at large r. Conversely, for < b < 1, the dissipation rate tends to the 
injection rate very slowly, oc t~^~ b \ and the energy is thus increasing with time as t b . 
The structure function S^ir) grows now at large distances as the zero mode r^ 2 ~^ b . For 
b > 1, coinciding particles do not separate and, in fact, separated particles collapse in 
a finite time. The consequences are that the dissipative anomaly is absent and that 
the energy is entirely transferred toward larger and larger scales in an inverse cascade 
process. The threshold 6=1 corresponds to the crossing of the exponents: (2 — £)6 
of the non-constant zero mode and 2 — £ of the constant-flux-of-energy solution. The 
picture is the mirror image of the one for the direct cascade, with the first exponent 
controlling now the small scale behavior and the second one appearing at the large 
scales. A sketch of the three possible situations is presented in Fig. 2. 

Concerning higher order correlations in the strongly compressible phase 6 > 1, we 
have shown that the inverse energy cascade is self-similar, i.e. without intermittency. 
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The effects of a large scale friction reintroduce, however, an anomalous scaling of the 
structure functions that do not thermalize without friction. They were exhibited in the 
explicit expressions for the p.d.f.'s of the tracer differences. As for the scalar density, 
the different behaviors of the Lagrangian trajectories were shown to result in the inverse 
or in the direct cascade of the total mass squared in, respectively, the weakly and the 
strongly compressible phase. As explained in Introduction, we expect the explosive sep- 
aration of the Lagrangian trajectories and/or their collapse to persist in more realistic 
ensembles of fully turbulent velocities and to play a crucial role in the determination of 
statistical properties of the flows at high Reynolds numbers. 

Acknowledgements. M. V. was partially supported by C.N.R.S. GdR "Mecanique 
des Fluides Geophysiques et Astrophysiques" . 



Appendix A 

Let us consider the operator M 2 of Eq. ( |2.16| ) on the half-line [r , oo[, r > 0, with 
the Neumann boundary condition d r f(r ) = 0. By the relation ( |2.1S| ), this means that 
we have to consider the operator N 2 on [uq, oof with the boundary condition 

d u u b -i v(u)\ a=UQ = (A.l) 

for u = r 2 . For non-integer b, the corresponding eigen-functions of N 2 are 

<Pe,u {u) = Ci(u ) J- h {u) + C 2 (u ) J b (u) (A.2) 



with u = yE/Z'u, see Eq. ( |2.20|) . Since 

M*) = Fr^)^( 1 + ^ 2 ))' ( A - 3 ) 
the boundary condition ( |A.1| ) implies that 

Ci(«o)0(uo) + C 2 (uv)0(uf~ l ) = 0. (A.4) 
As a result, for (non- integer) b < 1, lim S 2 {"°? = so that in the limit we obtain the 



eigen-functions of the operator N 2 . For (non- integer) b > 1, however, lim c ^[„°j = 

and the eigen-functions tend to those of N 2 + . The extension to the case of integer b is 
equally easy. 



Appendix B 

Let us give here the explicit form of the integrals of the kernels e~ tM z (r, r') against 
powers of r' . A direct calculation shows that for /i > and, respectively, b < 1 and 
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b > -1 

oo 



r(i-6) 



i^C-aSe' 1 - 6 '"^' (R1) 



• 1 ^ 1 (-^ + M + &,-§^). (B.2) 

A direct calculation gives also the asymptotic expansion of the kernels e~ tM ^(r, r') at 
small r: 

( t-d+2-t 
e 2 (r, r ) = J 



^iI(lT6)-(jT6)(4Z't)i+J (d*)J 1 - 



3=0 

An expansion around r' = may be obtained similarly. 



4Z< t 



(B.3) 



Appendix C 

Here we shall consider the long time behavior of the integral of the heat kernel of 
the operator M^: 

X(t,r,r') = fe- sM 2{r,r')ds. (C.l) 
o 

Using the explicit form ( p.22j ) of the heat kernel e~ sM 2 (r, r'), we may rewrite the last 
definition as 

t oo 

X{t,r,r') = fdsfe- sE E- b G(E;r,r')dE, (C.2) 
o o 

where 

G(E,r,r') = ^ ^ E b J^^E/Z'r^) J^^E/Z't'^) r'~ d+ ^ ~*. (C.3) 

Note that, by virtue of the relation flO|) , G(0,r,r') = r'~ d+1+a ~ e = G (r') 

and is independent of r. For finite times, the integration by parts, accompanier by the 
changes of variable sE <-> E gives the following identity: 

oo t oo 

bX(t, r, r') = t b f e- E E- b G(r 1 E, r, r')dE + J ds f e~ sE E 1 ^ 'd E G(E, r, r')dE. (C.4) 

0" 
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For b < 0, the t — > oo limit of X(t, r, r') exists and defines the kernel (M 2 



— \— l, 



r, r 



lim bX(t,r,r') 

t — >oo 



oo 




-6-1 



G(E,r, r') dE 

( i-d+2-t, 



for r < r' , 
for r > r' , 



(C.5) 



compare to Eq. (|3.4|) . The last expression has been obtained by the direct integration 
and the condition b < was required by the convergence at zero of the i?-integral (the 
kernels (M^) _1 (r, r') may be also found easily by gluing the zero modes of M 2 ). Note 
that the right hand side is a real analytic function of b = |^|. Now Eq. flC.4|) implies 
that 



lim 

t— >oo 



bX{t,r,r') - T(l-b)G (r')t b 



J E- b d E G{E,r,r')dE (C.6) 
o 



exists for b < 1 and is also real analytic in 6. On the other hand, it coincides with 
the long time limit in Eq. (|C.5 ) for b < and, consequently, must be given by the 
right hand side of this equation for all b < 1. The same arguments work after the 
integration of the above expressions against x{ r ')dpd{f') and hence the convergence of 
the expression ( |3.7| ) when t — > oo. 



Appendix D 

It is easy to give the exact expressions for the stationary 2-point functions of the 
scalars 6 and p for p < d ~ 2 2 g +e in the presence of positive diffusivity k. They are 



F 8 2 (r) = ^j-Jf^g^x^dUr') + ^ g(r) J f(r') X (r') dfi d (r') , (D.l) 



F 2 P (r) = f(r) J g(r') X (r') dp d (r') + ^1— /(r) g(r) / x(r') ^(r') , (D.2) 



where 



-d+l+q-g oo 

/(r) = (r« + f) ; and ff (r) = / /(C)" 1 £^ dC- (D.3) 



In the limit k — > these expressions pass into Eq. (|3.4j ) and ( |7.4|) , respectively. For 
d ~ 2 2 ; +s < p < similarly as at k = 0, the 2-point functions F 2 (t, r) and F 2 (t,r) are 
pumped into the constant and the f(r) zero modes of M 2 and (M 2 )*, respectively, and 



do not reach stationary limits, although the 2-point structure function of the tracer 
A 

e 



does. For p > j% the pumping into the zero modes reaches a constant rate. In the limit 



k — > 0, the zero mode f(r) into which F 2 is pumped goes to r ~ d+1+a -t for p < ^ but 
becomes the delta-function 5(r) for p > p. 



30 



Appendix E 



Let us prove the explicit expression ( |5.3|) for the even structure functions in the 
strongly compressible regime. Note that 

S e N (t,r) = £ (-l) IQC| ^(t,(rU,(0) 9eQC ) (E.l) 

QC{1,...,N} 

where Q c stands for the complement of Q. In the strongly compressible phase, by the 
multiple application of Eq. ( |5.1|) , we infer that 

^• S (W 9eQ) (0) qeQC ;r:) = J Pt s (v,0;r'(,v^ 

■ n - o n w - o dr " d A- (e.2) 

Substituting Eq. (|4.7|) into the expression (|E.1| ) and using the relations ( E.2|) , we obtain 

t 

S 9 N (t,r)= £ [ds [p£ s (r,0;T»,r'£ 

l<n<m<N q 

• ( E (-i) m K-M (ri) 96QX{n , m} , (^) ?eQC ) x(0) 

V {n,m}cQC{l,..,JV} 

+ E (-1) IQ1 (r'i'U X{ „ } , (rS), 60 c N{m} ) XlK - r£|) 

neQc{l,...,iV}\{m} 

+ E ("I) 101 WW,, K'W XM ) X(W - r' 2 '|) 

m6Qc{l,...,AT}\{n} 

+ E (-1) IQ1 W'U, (r 2 % eQCXM ) X(0)) drfdri' 

Qc{l,...,A r }\{ri,m} 7 

t <X> . /1 

= iV(iV-l) fdsjP*' s (r,r') S e N 2 (s,r') ( X (0)- X (r')) d^(0, (E.3) 


which is the sought relation. 
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